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Abstract 

A cross-layer optimization approach is adopted for the design of symmetric random access wireless 
systems. Instead of the traditional collision model, a more realistic physical layer model is considered. 
Based on this model, an Incremental Redundancy Automatic Repeat reQuest (IR-ARQ) scheme, tailored 
to jointly combat the effects of collisions, multi-path fading, and additive noise, is developed. The 
Diversity-Multiplexing-Delay tradeoff (DMDT) of the proposed scheme is analyzed for fully-loaded 
queues, and compared with that of Gallager tree algorithm for collision resolution and the network- 
assisted diversity multiple access (NDMA) protocol of Tsatsanis et at. The fully-loaded queue model 
is then replaced by one with random arrivals, under which these protocols are compared in terms of the 
stability region, average delay and diversity gain. Overall, our analytical and numerical results establish 
the superiority of the proposed IR-ARQ scheme and reveal some important insights. For example, it 
turns out that the performance is optimized, for a given total throughput, by maximizing the probability 
that a certain user sends a new packet and minimizing the transmission rate employed by each user. 

I. Background 

We consider a random access system with symmetric users who compete to communicate with 
a common receiver, or a base station (BS). Traditional approaches for analyzing such systems 
use the simplified collision model ([15] and references therein), which assumes that a message 
is received error-free by the BS if and only if a single user transmits. Under this model, several 
protocols, which attempt to avoid collisions, have been proposed in the literature, for example, 
Gallager tree algorithm (GTA) [8]. The collision model, however, does not adequately capture 
some important characteristics of the wireless channel, e.g., multi-path fading, and ignores 
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certain physical layer (PHY) properties like multi-packet reception (MPR) [10]. Recently, several 
researchers have started to focus on cross-layer optimization approaches which leverage the 
wireless medium to improve the performance of random access systems. For example, Naware 
et al. [10] analyzed the stability and average delay of slotted- ALOHA based random access 
channels with MPR at the BS. This analysis, however, abstracts out the physical layer parameters 
by using a very simplified model for MPR probabilities. Another example is [2] where Tsatsanis 
et al. propose the Network-assisted Diversity Multiple Access (NDMA) protocol, which uses a 
repetition based Automatic Repeat reQuest (ARQ) approach for collision resolution. As argued 
in the sequel, this protocol suffers from a significant loss in throughput resulting from repetition 
coding. In [1], Caire et al. analyzed the throughput of incremental redundancy (IR)-ARQ for 
the Gaussian collision channel with fully-loaded 1 queues and single-user decoders at the base 
station. By adopting the fully-loaded queuing model, this work ignores the stability issues that 
arise in practical random access systems with random arrivals. Moreover, the single-user decoders 
used in this work are sub-optimal and result in considerable throughput losses. To overcome the 
limitations of these previous works, we adopt a more realistic model for the physical layer, and 
develop a variant of IR-ARQ protocols optimized for random access over wireless channels. 

II. ARQ Random Access 

In this section, we introduce our system model and briefly review two existing random access 
schemes; namely, the GTA and NDMA protocols. To the best of the authors' knowledge, these 
two approaches represent the state of the art in the design of random access protocols. We then 
present an IR-ARQ random access protocol that overcomes the limitations of these protocols. 

A. System Model 

We consider a if-user symmetric random access channel with M antennas at each user and 
TV antennas at the receiver (base station). We assume that all the channels are independent and 
experience Rayleigh-flat and long-term static block fading where the channel coefficients remain 
constant during one collision resolution (CR) epoch and change independently from one epoch 

'Each queue has infinite packets for transmission. 
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to the next (a CR epoch will be defined rigorously in the next section). The channel coefficients 
are assumed to be perfectly known to the BS, but unknown to the random access users. We 
consider individual power constraints on the users, and denote the average received signal-to- 
noise ratio (SNR) of each user by p. In our model, time is slotted and a slot is composed of T 
channel uses. In order to control the number of users colliding in any particular slot, each user 
selects a slot for transmitting a new packet according to the probability-^ rule: in every slot, 
each user with a non-empty queue transmits a packet with probability p t and does not transmit 
with probability 1 — p t , where < p t < 1. We assume that the BS can perfectly identify the 
set of active users (by assigning a different control channel to each user). We initially assume 
fully-loaded queues in Section [Oil and then relax this assumption and consider a queuing system 
with random arrivals in Section |IVJ 

B. Gallager Tree Algorithm (GTA) 

This algorithm was proposed by Gallager [8] for the random access channel under the sim- 
plified collision model. The extension of this algorithm to our channel model mainly includes 
the probability -p t rule and an explicit assumption that the base station does not try to decode 
in the case of a collision. We describe the extended GTA as follows. The traffic in the channel 
is interpreted as a flow of collision resolution (CR) epochs. At the beginning of a CR epoch, 
each user uses the probability-^ rule to decide whether it should (or should not) transmit in that 
epoch. If none of the users choose to transmit, the slot remains idle and a new CR epoch starts 
from the following slot. If only one user chooses to transmit, then the message is assumed to be 
successfully decoded at the BS, and a new CR epoch begins from the following slot. But when 
a collision occurs, i.e., more than one user chooses to transmit in the current slot, the system 
enters into a CR mode, and only the users that participated in the collision at the beginning 
of a CR epoch are allowed to transmit until the end of that CR epoch. The colliding users are 
randomly split into two different groups according to a fair random split, wherein each user has 
an equal probability of joining either of the groups. A CR epoch is finished when all the users 
who have initiated it and not been excluded (or pruned) by the tree algorithm, obtain a slot to 
transmit their packets without collisions. We omit the detailed description of the algorithm for 
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brevity and refer interested readers to [8], [9]. 

C. Orthogonal Network-Assisted Diversity Multiple Access ( O-NDMA ) 

The NDMA protocol was proposed by Tsatsanis et al. [2] and relies on the use of time 
diversity through a repetition ARQ scheme to resolve collisions between users. At the beginning 
of each CR epoch, the transmission of each user will be determined by the probability-p t rule 
as in the GTA protocol. If none or only a single user choose to transmit, then the next CR 
epoch starts from the following slot as before. However, when k (> 2) users transmit, then all 
those users repeat their transmissions in the next (k — 1) slots. At the end of k slots, the BS 
is assumed to be able to reliably decode the k packets, and a new CR epoch begins from the 
next slot. On the other hand, in [3], Zhang et al. proposed a new variant of NDMA which does 
not rely on time diversity to resolve/detect collisions. This variant, named B-NDMA, relies on a 
blind signal separation method utilizing a Vandermonde mixing matrix constructed via specially 
designed user retransmissions. In B-NDMA, the detection and resolution of a A;-user collision 
require (k + 1) slots. However, in this paper, we assume the use of separate control channels for 
collision detection; which allows for a slightly more efficient version of the B-NDMA protocol, 
named orthogonal NDMA (O-NDMA), which requires only k slots to resolve a /c-user collision, 
without relying on temporal diversity. The behavior of users in O-NDMA is the same as that in 
NDMA, with the only difference that in case of a /c-user collision, user i transmits its symbols 
scaled by (wi) e = (e^Y, where i = 1, • • • , k and j = \f—l, in the £-th slot after the initial 
collision. At the end of the k th slot, the BS utilizes the orthogonal structure constructed by w^s 
to decompose the joint decoding problem into k single-user problems. For example, suppose 
that user 1 and user 2 have collided (k = 2), and user i's codeword is x i5 for i = 1,2. Then, 
the BS coordinates the users so that user 1 repeats x x whereas user 2 transmits — x 2 , in the slot 
following the collision. To decode user 1, the BS calculates the sum of the received vectors in 
the two slots, while to decode user 2, it takes the difference (i.e., matched filtering). This way, 
the multi-user interference is removed, and single-user decoders can be utilized to recover both 
packets. It is worth noting that O-NDMA requires symbol-level synchronization to facilitate the 
interference cancellation described above. Hence, our results for O-NDMA can be interpreted 
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as optimistic upper bounds on the performance of repetition based random access protocols. 

However, O-NDMA is still sub-optimal for two reasons. First, the BS might be able to decode 2 
the messages of k colliding users in less than k time slots. Conversely, it is also possible that 
k time slots are insufficient for the successful decoding of the k packets. Thus, such a static 
strategy may result in a throughput loss. Second, O-NDMA is essentially a repetition based 
collision resolution mechanism. Although this results in a low-complexity decoder at the BS, the 
throughput performance is highly sub-optimal, as shown rigorously in the sequel. A significant 
improvement in the throughput can be achieved by allowing for IR transmissions from the 
colliding users within the CR epoch, and using joint decoding, across ARQ rounds and users, 
at the base-station (as discussed next). 

D. IR-ARQ Random Access 

To overcome the disadvantages of the existing protocols, we propose a new IR-ARQ random 
access protocol operating as follows. Each user encodes an information message (packet) of 
B T bits using a codebook of length-LT codewords, where L is an integer denoting a deadline 
constraint on the transmission delay (i.e., a constraint on the maximum number of allowed ARQ 
rounds). Codewords are divided into L sub-blocks of length T. At the beginning of each CR 
epoch, the users choose to transmit or not based on the probability-^ rule as before. Once a 
user chooses to transmit in a particular slot, it transmits its first T symbols during that slot. 
On receiving signals at the end of a slot, the BS uses a joint decoder that decodes the received 
observations both across users and ARQ rounds. If the receiver successfully decodes all the 
transmitted messages, it feeds back an ACK bit; otherwise, it returns a NACK signal. On receiving 
an ACK, the CR epoch is terminated and a new CR epoch starts from the next slot. Thus a CR 
epoch can be defined as the time between two successive ACKs from the receiver (we observe 
that this definition requires the BS to return an ACK message after an idle slot). On the other 
hand, if a NACK is received, each colliding user sends its second sub-block of T codeword 
symbols in the next slot, while all the other users remain silent. The ACK/NACK rule applies in 

2 Multiple messages can be jointly decoded in a single transmission block, with an arbitrary small error probability, if a 
rate-tuple lies within the capacity region of the channel and a sufficiently large block length is used [14]. 
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a similar manner, until the L th slot is reached (after (L — 1) consecutive NACKs). In this case, 
the receiver sends an ACK regardless of its decoding result. If a certain message is decoded 
after £ ARQ rounds, the effective coding rate for the corresponding user becomes R/£ bits 
per channel use (BPCU), where R = (Bt/T) denotes the rate computed assuming only one 
transmission round. Finally, we note that, unlike the O-NDMA, the IR-ARQ protocol requires 
only slot-synchronization. 

III. The Diversity-Multiplexing-Delay Tradeoff (DMDT) 

In this section, we analyze the DMDT of the proposed IR-ARQ protocol and contrast it with 
our two benchmark protocols under the assumption of fully-loaded queues. The "fully-loaded" 
assumption allows for analyzing the maximum achievable throughput without focusing on the 
stability and average delay issues, for the moment. 

A. Definitions 

We borrow the notion of DMDT from [5]. This notion is a generalization of the Zheng- 
Tse diversity-multiplexing tradeoff (DMT) which characterizes the fundamental limits of fading 
channels in the high SNR regime [6]. The delay here refers to the maximum transmission 
delay corresponding to our upper bound L on the number of ARQ rounds (including the first 
one). In particular, we consider a family of ARQ protocols where the size of the information 
messages B T (p) depends on the operating SNR p. These protocols are based on a family of 
space time-codes {C p } with a first round rate of R(p) = B T (p)/T and an overall block length 
TL. For this family of protocols, we define the first round multiplexing gain r and the effective 
ARQ multiplexing gain r e as 

1- R {P) A A y VFL(P) n , 

r = hm — —- and r e = hm — . (1) 

Here 7] F l(p) is the average throughput of the ARQ protocol in the random access channel with 
Fully-Loaded (FL) queues, i.e., 

VFL (p) = hm (2) 
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where s is the slot index and b(s) is the total number of message bits transmitted up to slot s. 
Note that the message bits received in error at the BS are also counted in b(s). The effective 
ARQ diversity gain is defined as 

d = - lim lQ f Pe(p) , (3) 

p^oo log 2 p 

where P e (p) is the system error probability, which is defined as the probability that at least one 
of the messages is not correctly decoded by the BS. In the symmetric random access channel, 
the diversity gain obtained from © is the same as the diversity gain of an individual user, since 

K 

P eW (p)<P e (p)<5> ea) (p), Vie{l,---,K} , (4) 

[7] where P e (i) (p) is the error probability of the i th user. In summary, the DMDT of a certain 
protocol characterizes the set of achievable tuples (d, r e , L) (here, we observe that our results are 
information theoretic in the sense that we assume the use of random Gaussian codebooks [6]). 

In our analysis, we will make use of the results of Tse et al. on the diversity-multiplexing 
tradeoff of coordinated multiple access channels [7], where the access mechanism is controlled 
by the base-station. In the sequel, we denote the diversity gain of the coordinated multiple access 
channel with k users as d^ AC (r), which is given by 

( d M ' N (r), r<min{M, I ^ T } 
d™ AC {r) = { ^ h ~ 1 ,k+li , (5) 

[ d kM > N {kr), r > min{M, ^} 

where d M ' N (r) is the diversity gain of the point-to-point channel with M transmit and N receive 
antennas, and multiplexing gain r, as given in [6]. 

In the ARQ setting, we denote the event that a NACK is transmitted in the £ th ARQ round, 
when k users are transmitting simultaneously, by Ak{£), for £ = 1, • • • ,L — 1, and the error 
event in the L th round by Ak{L). We also denote the complement of Ak(£) by Ak{£)- We 
define a k (£) 4 Pr (^ fc (l), ■ ■ • , A k (£ - l),A k (£)) and (3 k (£) 4 Pr (A k (l), ■ ■ • , A k (£)) for 
£ = 1, ■ ■ ■ , L, where, by definition, we let (3 k (0) = 1, for k = 1, • • • , K. Note that a k {£) is the 
probability that the length of a CR epoch is £ (slots), given that k users have collided initially. 
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Following the approach of [2], we classify the CR epochs from the viewpoint of a particular 
user (say user 1) into either relevant or irrelevant epochs, depending on whether a packet of that 
particular user is being transmitted in that CR epoch or not. The lengths of the relevant and the 
irrelevant epochs of user 1 are random variables, which are denoted by U and V, respectively. 
For notational convenience, we denote the pmf of a Bernoulli random variable with population 
K and probability of success p by, B(K, k,p) = (, ) p h (l — p) K ~ k . 

B. Main Results 

First, we characterize the DMT for GTA (note that we do not have a deadline in this protocol, 
and hence, no limit on the maximum transmission delay). 
Proposition 1: The DMT for GTA with a given p t e (0, 1] is 



where X k and J k can be found by the following recursions: 

k 

X k = 1 + B(k, 0, 0.5)X k + B(k, 1, 0.5) (1 + X k - X ) + B(k, i, 0.5)^ , (7) 

i=2 

k 

J k = B(fc,0,0.5)J fe + B(fc,l,0.5)(l + + ^2B(k,i,0.5)Ji , (8) 

i=2 

for k = 2, 3, • • • , with X = X\ = 1 and J = 0, J\ — 1. 

Proof: Noticing that the CR epoch termination event is a renewal event under the fully- 
loaded assumption, the result can be easily derived by extending the recursion analysis in [9] 
and using the renewal-reward theorem [16]. The details are omitted due to space limitation. ■ 
Since the GTA protocol is inspired by the simplified collision model, the main idea is to assign 
a single slot exclusively for transmission of each colliding user (that was not pruned by the 
algorithm). The resulting DMT, therefore, is given in terms of a single-user performance, i.e., 
d\ IAC {.). The main drawback of the algorithm is the relatively large number of slots needed to 
resolve each collision, which translates into a loss in the effective multiplexing gain, i.e., the 
argument of d^ AC (.) in ©. It is now easy to see that GTA cannot achieve the full effective 
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multiplexing gain of the multiple access channel, i.e., mm{KM, N}. An example highlighting 
this fact will be provided in the later part of this section. On the other hand, the DMT in © 
reveals the performance dependence on p t (and r), which implies the possibility of maximizing 
the diversity gain by choosing the appropriate values, p* t and r* for each r e 6 [0, mmjfTM, N}). 
At the moment, we do not have a general analytical solution for this optimization problem. 
However, the solution for the special case of two users is obtained in Section IIII-CI 

Next, we characterize the optimal DMT for the O-NDMA protocol (Again we do not have a 
delay parameter in the tradeoff since the number of ARQ rounds is always equal to the number 
of colliding users). 

Proposition 2: The optimal DMT for O-NDMA is, 

d ONDMA (r e ) = df IAC (r e ). (9) 
Proof: The DMT for O-NDMA with a given p t G (0, 1] and r is found as 

dONDMA^ = d MAC( r) f = K Pt + (1- Pt ) K ^ 

utilizing the average throughput results in [2], and noting that the average SNR of each single- 
user decoder is still p. Then, it is easy to find that the optimal values (r*,p^) = (r e , 1), which 
yields ©. We omit the details due to space limitation. ■ 
The matched-filter-like structure utilizing the orthogonality of transmissions over different slots 
allows the O-NDMA protocol to achieve the single-user performance, as we see from (ITOT) . 
Furthermore, p* t = 1 ensures that the throughput is maximized, and the optimal DMT is given by 
©. By comparing the expressions in © and ©, we realize that the O-NDMA protocol achieves 
a larger diversity gain, as compared with the GTA protocol, for any r e less than min{M, iV}. 

Finally, the optimal DMDT of the IR-ARQ random access protocol is characterized in the 
following theorem. 

Theorem 3: The optimal DMDT for the IR-ARQ protocol is, 

d IR (r e ,L) = < fAC (^). (ID 
Proof: (sketch) First, we assume an asymptotically large block length T — >• oo to allow our 
error correction (and detection) scheme to operate arbitrarily close to the channel fundamental 



limits. An application of the renewal-reward theorem [16] gives 

PtKR no , 

Vfl = t? rn • (12) 

1 + Ef=i B ( K > k, Pt) Eti Phi® 

In addition, given that joint typical-set decoders [1], which have an inherent ability to detect 

errors, are used for the channel output over slots 1 to £ in ARQ round I, extending the results 

in [5] and [7], the probability of error P e is upper-bounded by, 

K 

P e < Y, B i K ^Pt)^i L )- (13) 

fc=l 

Noting that in the high SNR regime f3 k (£) approaches 

( ( £N}\ a I 0, r<min{£M, ^} 

\im p k (p,£) = l(r>mml£M,— \) ± { . (14) 



h 1 i 1 1, r > min{£M, f^} 



we find the DMDT with a given p t e (0, 1] as 



d IR (r e ,L) = d% AC (^\, (15) 



where r can be obtained from r e using the relation (for < r < min{M, N}\ 

p t Kr 



(16) 



1 + Eti B(K, k,p t ) Ef=i 1 (r > min{/M, f }) 
from ([T2HI1I>, and the results in [5], [7]. Finally, we find the optimal values [r*,p* t ) = 1), 
which gives dTTT) . A detailed proof is provided in Appendix U ■ 
Two remarks are now in order. First, we elaborate on the intuitive justification for the optimal 
values (r*,pl) = (^,1) for the IR-ARQ protocol. In the asymptotic case with p — > oo, the 
error probability is dominated by the worst case f^-user collision for any p t E (0, 1], which does 
not depend on p by definition. This implies that choosing p t = 1 will maximize the average 
throughput, without penalizing the asymptotic behavior of the error probability. Now with p t = 1, 
choosing r* = ^ < min{M, ^} will result in an effective multiplexing gain equal to r e and 
will minimize the number of rounds needed to decode the colliding messages, since each user 
is transmitting at a small rate. Furthermore, it is clear that with this choice of r* , we can 
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achieve any desired effective multiplexing gain less than min{KM, N} (the degrees of freedom 
in the coordinated multiple access channel). Next, Comparing Propositions [T] El and Theorem 
|3J it is straightforward to verify that the DMT of the IR-ARQ protocol is always superior to 
that of the GTA and O-NDMA protocols. This advantage of IR-ARQ is a manifestation of 
the ARQ diversity resulting from the IR transmission and joint decoding. More specifically, 
the ARQ diversity scales down the effective multiplexing in the right hand side of dTTb . and 
hence, results in an increased diversity advantage (since dj^ AC (.) is a decreasing function in its 
argument). The O-NDMA protocol does not allow for efficiently exploiting the ARQ diversity 
due to the sub-optimality of repetition based ARQ. 

C. Examples 

We numerically illustrate the gains offered by the IR-ARQ protocol, as compared with the 
GTA and O-NDMA protocols, in the following two examples. 

1 ) Two-User Scalar Random Access Channels: We consider the single- antenna 2-user random 
access channel, i.e., M = N = 1 and K = 2. Substituting these parameters in Proposition [H we 
obtain the DMT for the GTA protocol as, d GTA (r e ) = df AC (^^r}j = 1 - (^Pj r e . In 
order to maximize the effective multiplexing gain that achieves nonzero diversity gain, we need to 
choose p t = 4g, which yields the optimal DMT for GTA as d GTA (r e ) = 1 — \/3r e , < r e < 
The optimal DMTs for O-NDMA and IR-ARQ are obtained from Proposition |2] and Theorem 
13 Fig. [T] compares the tradeoffs of the three protocols where the IR-ARQ protocol is shown to 
dominate our two benchmarks, with both L = 1,2. Even though O-NDMA achieves the nominal 
single-user DMT without multi-user interference, i.e., d(r e ) = 1 — r e , Vr e < 1, it is still worse 
than IR-ARQ, since it wastes slots to facilitate single-user decoding and relies on the repetition 
ARQ. In addition, as L increases from 1 to 2, the DMDT of IR-ARQ improves, as expected. 

2) Two-User Vector Random Access Channels: We consider a 2-user vector random access 
channel with M = 1 and iV = 2. By allowing multiple antennas at the BS, the total degrees of 
freedom of the system is increased by a factor of 2, as compared with the scalar channel in the 
previous example. The tradeoffs achieved by the three protocols in this scenario are shown in 
Fig. El First, we observe that the three protocols achieve an increased diversity gain, for a given 
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r e , when compared with the scalar channel in Fig. [T] However, the full effective multiplexing 
gain, r e = 2, is not achieved by the GTA and O-NDMA protocols, since these two protocols 
exclude the possibility of first-round decoding when a collision occurs. The IR-ARQ protocol, 
on the other hand, achieves r e = 2, and the DMDT further improves as L increases. 



In this section, we relax the "fully loaded" assumption adopted in Section |TTT| In addition to the 
traditional measures of stability region and average delay commonly used in this set-up, we also 
consider the probability of error. In particular, for the proposed IR-ARQ protocol, we will show, 
through numerical results, that the choice of the transmission-delay constraint L determines an 
interesting tradeoff between the average delay and error probability: For typical SNR, a larger 
L leads to an increase in the average delay along with a decrease in the error probability. 

We consider infinite-length queues at the users, that are fed by randomly-arriving packets 
of a fixed length of Ba information bits. For simplicity, we assume that Bt = Ba = B, 
i.e., the arrival packet size and the transmission packet size are the same. Thus the first-round 
transmission rate R is equal to the arrival rate Ra = (B/T). To emphasize that R is a system 
parameter determined by the arrival packet size, we denote the first-round multiplexing gain r 
by ta in this section, and call it the arrival multiplexing gain. The packet arrival rate of user 
i is Aj = X/K packets/slot, where A denotes the total packet arrival rate, where arrivals are 
assumed to be independent across users. 

A. Stability and Average Delay 

We use the following notion of stability [4]: let g(m) = (gi(m), ■ ■ ■ , <?k-(to)) be the vector 
of the backlogs at the beginning of CR epoch m. Then, queue i of the system is stable if 
linim^oo Pr(<7j(m) <g) = F(g) and ]im.g^ 00 F(g) = 1. Furthermore, we say that the system 
is stable if all the K queues in the system are stable. 

The stability region of GTA and O-NDMA can be found using the techniques in [9], [4] as 



IV. Random Arrivals 



A < 



and 



A < 



K Pt + (1 - p t ) 



Kp t 



(17) 
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From the literature ([10] and references therein), we find only limited analytical results on the 
average delay of slotted ALOHA channels. In this paper, we present only numerical results 
for the average delay of the GTA and O-NDMA protocols, and provide an approximate delay 
analysis for the proposed IR-ARQ protocol. The average delay of the IR-ARQ scheme can 
be approximated by using the analysis of the M/G/l queue with vacations [15], following the 
approach of Tsatsanis et al. [2]. This analysis yields only an approximation of the average 
delay, since the CR epoch lengths of the IR-ARQ scheme are dependent on the traffic load (and 
hence are not independent and identically distributed (i.i.d.) as needed for the result to hold). 
However, as we will see, the i.i.d. property holds in the limit of p — > oo, and hence, our result 
becomes asymptotically accurate. We also note that as K increases, this approximation becomes 
progressively more accurate [2]. We summarize our results in the following theorem. 

Theorem 4: Assuming that 3 i < oo with £ E {1, • • • , L}, such that olk[£) > 0, the neces- 
sary and sufficient condition for the stability of the IR-ARQ protocol is (in packets/slot) 



Vfl PtK 



For Poisson arrivals, when A satisfies (fT8t . the average delay is approximately given by (in slots) 



/I \ A (E t/ 2 + (2 - p f- pt W 2 ] + 2 - l) E[UE[v) m V 2] 
D « E[U}+( - - 1 E[V] + ^ 7 7 ^ 4- ~ + , 

U J 2(k-\(e[u] + (±-i)e[v})) 2E[y]' 

(19) 

where the expected values are evaluated as, 

K L-1 K L-1 

e[u) = i + J2 B ( K - x > k - X 'P) J2 p*w ; E [ [/2 ] = 1 + 12 B{ - K - 1; k - 1 ' p) 12 (2i + 

k=l i=l k=l e=i 

K-l L-1 K-l L-1 

nV} = l+Y l B(K-l,k,p)J20kW ; E[y 2 ] = l+^ 6(^-1, fc,p)^(2£+l) / 3 fe (^), 

k=l £=1 fe=l 1=1 



and p e (0, 1] satisfies Kp = A [l + £f =1 B(K, k, p) Zti PkW 



(20) 
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Moreover, the delay expression in (fT9l) holds with probability 1 if U and V are i.i.d. and U and 
V are mutually independent. 

Proof: See Appendix |TT| ■ 
A few remarks are now in order: First, the assumption in Theorem |4] always holds when L is 
finite since the length of any CR epoch is bounded by L. If L = oo, this assumption requires 
the existence of a non-zero probability that the length of an epoch is finite. Second, as p — > oo, 
the stability region (fTSt approaches 

x < VjF_ 

1 + Ef=i Eti 1 1 (r A > min {£M, f }) ' 

To achieve the maximum stability region, we need to maximize the right hand side of <|2TJ over 

p t . At the moment, we do not have a general solution for this problem. Thus, we present results 

only for one special case: r A < min{M, ^}. In this case, the stability region is A < ptK, and 

the maximum stability region is thus given by A < K for the optimal choice of p t = 1. This 

is a remarkable improvement over the O-NDMA protocol, whose maximum stability region is 

only A < 1, for any r A - Finally, the diversity gain with random arrivals can be readily obtained 

from the results in the previous section. The only difference is that, unlike the fully-loaded case, 

one cannot optimize over r A in this scenario. In summary, we find that the GTA, O-NDMA and 

IR-ARQ protocols achieve the diversity gains d GTA (r A ) = d ONDMA (r A ) = df AC (r A ) and 

d IR (r A ) = d% AC (^), respectively. 

B. Examples 

1 ) Two-User Scalar Random Access Channels: Here, we consider the random access channels 
with M — N — 1. For ease of analysis, we assume that L > K = 2 for the IR-ARQ scheme. 

The stability region of the different random access protocols with p — ► oo is summarized in 
Table |U In addition, the error probability, diversity gain and average delay are shown in Fig. 
Fig. |4] and Fig. |5] respectively. Here, the stability region and diversity gain of the three protocols, 
and the average delay of the IR-ARQ protocol with p — ► oo, are obtained analytically. However, 
the average delay of the GTA and O-NDMA protocols, and the average delay of the IR-ARQ 
scheme with p < oo are obtained through numerical simulations. In these simulations, we use 
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R = 7~Alog(l + p) with r A = 0.45, and p t — 1 for the IR-ARQ and the O-NDMA protocols, 
while p t — 4g for the GTA protocol. It is assumed that the transmission results in errors, if and 
only if the channel is in outage [1]; which is a valid assumption if T is sufficiently large. In 
addition, for the IR-ARQ protocol, it is assumed that the errors in £ th round, where £ < L, are 
always detected. We also note that when r A < 0.5, the average delay expression for the IR-ARQ 
scheme, evaluated from Theorem |H holds with probability 1, and is given by (when p t = 1) 
D = 1.5 + 2 (^x) • Table U and Fig. |4| shows that both the stability region and diversity gain 
of the IR-ARQ protocol are the largest. Next, we focus on the delay and the error probability 
of IR-ARQ with different L's and different p's reported in Fig. |3] and Fig. |5j We observe that 
the delay approaches the asymptotic result with p = oo, and the difference of the delay for 
IR-ARQ with L = 2 and with L = 4 decreases, as p increases, which agrees with the analytical 
results. Furthermore, Fig. |3] and Fig. |5] reveal an important insight into the relation between the 
performance of IR-ARQ and the transmission-delay constraint L, i.e., a tradeoff between average 
delay and error probability emerges. These figures suggest that for certain finite p's, a large L 
achieves a small error probability, at the expense of a large average delay and a small stability 
region. Therefore, depending on quality-of-service (QoS) requirements, L can be adjusted for 
achieving the best performance. 

2) Two-User Vector Random Access Channels: Here, we consider the 2-user random access 
protocols with M = 1 and iV = 2 in the high SNR regime (p — > oo). We first note that 
the stability region and delay of the GTA and O-NDMA protocols are not different from the 
scalar case; only the diversity gain changes with this multiple-antenna setting. For the IR-ARQ 
protocol, on the other hand, the average delay is given by D = 1.5 + ^TaJ f° r an Y r A £ [0, 1), 
and the stability region is given by, A < 2, < r A < 1, with p t = 1. Comparing the stability 
region of the vector IR-ARQ protocol with that of the scalar IR-ARQ protocol, we find that the 
vector IR-ARQ achieves a better stability region, especially when r A > 0.5. Finally, Fig.[6]shows 
the diversity gain achieved with different random access protocols. As expected, the IR-ARQ 
protocol achieves the best diversity gain, which improves as L increases. 
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V. Conclusions 



We have proposed a new wireless random access protocol which jointly considers the effects 
of collisions, multi-path fading, and channel noise. The proposed protocol relies on incremental 
redundancy transmission and joint decoding to resolve collisions and combat multi-path fading. 
This approach represents a marked departure from traditional collision resolution algorithms and 
exhibits significant performance gains, as compared with two benchmarks corresponding to the 
state of the art in random access protocols; namely GTA and O-NDMA. It is interesting to 
observe that, in order to fully exploit the benefits of the proposed IR-ARQ protocol, all the 
users with non-empty queues must transmit with probability one, when given the opportunity, 
and should use a small transmission rate. Finally, we have identified the tradeoff between average 
delay and error probability exhibited by the IR-ARQ protocol for certain SNRs, and have shown 
that this tradeoff can be controlled by adjusting the maximum number of ARQ rounds. 

Appendix I 
Proof of Theorem 

To find the long-term average throughput of IR-ARQ, we first focus on the distribution and 
the expected value of the relevant and the irrelevant epochs for user 1, U and V. The probability 
mass functions (pmf) of U and V are 

K 

Pr(U = £) = ^B(K-l,k-l,p t )a k (£), (£=!,■■-, L), (22) 

k=l 



for any (a 1; ■ • • , a L ) 6 R L . Using this relation, it is straightforward to get 

K L K-l L 

E[U] = ^BfK-U-l.p^ftfM) ; E[V] = (l- Pt ) K - 1 + £ B(K-1, k, Pt ) ]T/3 fc (£-l). (25) 
fc=i 1=1 k=i e=i 
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(23) 



We introduce the relation shown in [5] for deriving the expected values of U and V: 




(24) 



Now, we are ready to calculate the long-term average throughput (IT2T> and the upper-bound of 
error probability P e given in ( fT3t for the IR-ARQ scheme as in the following. First, we prove 
([T2l utilizing the renewal theory [16]. Denoting the average throughput of user 1 by r]\ , the 
average throughput of the symmetric system is given by t}fl = Krji . Under the fully-loaded 
assumption, the event that a CR epoch terminates is a renewal event. We associate a random 
reward TZ to the occurrence of the renewal event; TZ — R BPCU if the CR epoch is a relevant 
epoch for user 1, and TZ — otherwise. Then, the renewal-reward theorem [16] with (l25t gives, 



.. h{s) E[K] p t -R+{l-p t )-0 Pt R 

//; = Inn 



sT E[X] Pt ■ E[U] + (1 - p t ) ■ E[V] 1 + £f =1 B(K, k, p t ) Pk{Q ' 

(26) 

Since T] FL = Kr\\, we obtain i] FL as given in (fL?! . Next, we prove (fT3t . An error occurs in the 
IR-ARQ system in two different cases: (i) when decoding failure is not detected at the BS and 
an ACK is fed back, or (ii) when decoding fails at round L. Let E k (£) denote the event that the 
decoder makes an error with i received blocks when k users have collided in the first round. 
Then, we can upper-bound P e as [5] 

K L 



= ^ B(K, k, p t ) Pr(E k (£), an epoch length when k users have collided = £) 

k=i e=i 

K K 

< Y< B ( K i k iPt) Pr(£fc(£),A,--- ,A L -i) + e = J2 B ( K ^ k ^t)Pk(L)+e, (27) 



k=l k=l 

where e — > as T — > oo. The intuition behind this upper-bound is that the undetected error 
probability approaches zero as T — > oo for the joint typical-set decoder, and hence the error 
probability is dominated by the information outage probability. 

On the other hand, the diversity gain <fT3T) can be found as in the following. We first find 
Pk(p, = p~ d k IAC ^ r / l \ using the results in [5] and [7], where A(p)=p~ b implies b = — Hindoo lo ^ j4( - / " 
as defined in [6] and <, > are similarly defined. With this, given that p t does not depend on p, 
(fT3t implies that P e (p) satisfies the exponential inequality P e (p) < 

p-< AC (f) 5 as t -»• oo. In ad- 

dition, applying the outage bound in [5] to the random-access system yields P e (p) > p K ^ T > ■ 
These two exponential inequalities imply (fT5t . 



\2P 
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Noticing that the span of d M ' N (r) is r G [0, min{M, N}) [6], we verify that the span of 
d% IAC (j) is r G [0,min{£M, ^jr}). Then (fT4T) can be readily verified as in the following. For 
r < min{£M, f^}, it is obvious that fi k (p,t) -> as p -> oo from /3 fc (p,f) = p- d k AC W). On 
the other hand, if r > min{£M, then (3 k (p,£) — > 1 as p — ► oo since the outage probability 
approaches 1 as p — > oo, and the error probability given the outage event approaches 1 as 
T — > oo by the strong converse [14]. Combining the results shown above, we get the DMDT of 
the proposed IR-ARQ protocol as (fT5l) . where the relation between r and r e is given in (fToT) . 

Finally, we consider the optimal pairs (r*,p£) that achieve the largest diversity gain for a 
desired r e . Regarding r e in (IT6b as a function of r, we observe that r e (r) is discontinuous at 
the points r = min{£M, ^}, £ = 1, • • • , L — 1 and k — 1, • • • , K. We consider the values 
of r that lie within the first discontinuity of r e (r), i.e., r G [0, min{M, ^}). For these r, ( fTBT ) 
yields r e = p t Kr. Since the diversity gain increases when r decreases, we want to determine 
the smallest value of r that achieves the desired r e . We find that r is minimized when p t = 1. 
Thus for r G [0, min{M, ^}), the optimal choice of (r,p t ) achieving r e is 1). Furthermore, 
we find that this choice achieves the maximum effective multiplexing gain given by the degrees 
of freedom of the channel (mmjfTM, N}). Thus we do not need to consider the values of 
r > min{M, ^}, since it is clear from (fToT) that such r values result in a smaller diversity gain 
for the same desired r e . 

Appendix II 
Proof of Theorem 

We consider the backlog evolution g(m) of IR-ARQ, where m is the epoch index. We observe 
that g(m) is an embedded Markov chain; gi(m), the backlog evolution of user i is, 

{(gi(m) — 1) + + cii(m), with probability p t 
(28) 
<7i(m) + aj(m), with probability 1 — p t 

where cij(m) is the number of packets that arrived at user i's queue during epoch m, and = x 
if x > 0, = otherwise, for a real number x. 

We first prove that (fT8t is the necessary and sufficient condition for the stability of IR-ARQ. 
One can straightforwardly prove that under the assumption that there is a finite £ with a K (£) > 0, 
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g(m) is a homogeneous, irreducible and aperiodic Markov chain, by following the argument in 
the proof of Proposition 1 in [13]. Given that the Markov chain has these properties, stability 
of the system is equivalent to the existence of a limiting distribution for the Markov chain, and 
thus is also equivalent to ergodicity of the Markov chain [11], [12]. Sufficiency and necessity of 
(TTSt for the ergodicity can be straightforwardly proved by following the footsteps of the proof 
of Theorem 1 in [11]. In particular, we consider a stochastically dominant system, in which 
users are first chosen according to the probability-p t rule, and those users with empty queues 
transmit dummy packets. It can be shown that (TTSt is a sufficient condition for the stability of 
the dominant system, which implies the stability of the original system. On the other hand, we 
observe that the bounded homogeneity property [12] holds for the Markov chain (l28t . and thus 
the instability of the dominant system implies the instability of the original system. 

Next, we consider the approximate average delay. We denote the time period between the 
instance when a packet of user 1 reaches the head of its queue and the instance when it has 
finally been transmitted, by a random variable y (slots). Then, the average delay D is given 
by the result for the M/G/l queue with vacations [15], D = E[y] + 2( - 1 AE ^^^ + ^SU, with 
probability 1, if y is i.i.d. and V is i.i.d.. Here, the first moment of 3^ is calculated as, 

E[y\ = E\p t U + Pt {l- Pt ){U + V)+ pt(l - Ptf{U + 2V) + ■■■}= E[U] + f i - E[V]. 

On the other hand, the second moment of y is calculated as, 

E[^ 2 ] = E[p t U 2 + p t (l - p t )(U + V) 2 + p t (l - Pt ) 2 (U + 2V) 2 + ■■■] 
= E[U 2 } + (2 - pt) J-}- pt) E[V 2 } + 2 (± - l) E[U]E[V], 

if U and V are independent. Substituting the values of E[3^ 2 ] and E[y] in D, the approximate 
delay ( fT9l can be readily obtained. The expected values of the steady-state epoch lengths, 
E[U], E[U 2 ], E[V] and E[V 2 ] used in the expression dT9l are easy to derive utilizing (l24b . noticing 
that the pmf's for U and V are given by (l2"2"b and (|2"31 with a substitution of p into p t , where 
p = p t (l — p e ) and p e is the steady-state probability of a user's queue being empty. 

Finally, to see that the steady-state transmission probability p can be found by solving the 
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equation d20l . we consider the method of the steady-state analysis of the Markov chain gi(m), 
whose time-evolution is given by d2"%t. Let us define the steady state values: g\ = linv^oo g\{m) 
and a x = linim^oo a x {m). Then, taking expectation on both sides of (|2"%1) results in: 

E[ gi (m + 1)] = E[ gi (m)} - p t Pr(^(m) > 0) + E[ai(m)]. (29) 

In the limit as m — ► oo, we have, E[#i] = E^] — PtPr(#i > 0) + E[ai], or p = E[oi]. Thus, 

p = ^(p t Pr(g 1 >0)E[U} + (l-p t Pt(g 1 >0))E[V}) = ^{pE[U} + (l-p)E[V}), (30) 

which is equivalent to (l20l) . 
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TABLE I 

Stability Region of Different Two-User Scalar random access Protocols 





Stability Region 


Maximum Stability Region 


GTA 


\<2p t /(l + 3p 2 t ) 


A < l/y/3, with Pt = l/y/3 


O-NDMA 


\<2 Pt /(2 Pt + (i- Pt y) 


A < 1, with p t = 1 


IR-ARQ 


\ 


\ A < 2p t , r A < 0.5, 
[ \<2p t /(l+p 2 t ), r A >0.5. 


1 


f A<2, r4<0.5, ., 1 
. ' _ _ , with p t = 1. 
A < 1, r A > 0.5. Ft 




0.3 ■ 
0.2 
0.1 ■ 



■GTA 
-O-NDMA 
-V— IR-ARQ, L=1 
IR-ARQ, L=2 
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Effective Multiplexing Gain r 



1. Diversity-multiplexing tradeoff for various two-user scalar random access systems. 
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Fig. 2. Diversity-multiplexing tradeoff for various two-user vector random access systems. 
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Fig. 3. System error probability versus SNR for various two-user scalar random access systems with random 
arrivals. Here, ta = 0.45. 
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Fig. 4. Diversity gain versus the arrival multiplexing gain Ta for various two-user scalar random access systems 
with random arrivals. 
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Fig. 5. Load per user versus average delay for various two-user scalar random access systems with Poisson arrivals. 
Here, r A = 0.45. 
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Fig. 6. Diversity gain versus the arrival multiplexing gain r& for various two-user vector random access systems 
with random arrivals. 
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